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ABSTRACT 

We present new XMM-Newton observation of the Ultra Luminous Infrared Galaxy (ULIRG) NGC 6240. 
We analyze the reflecting grating spectrometer (RGS) data, and data from the other instruments, and find a 
starburst dominated 0.5-3 keV spectrum with global properties resembling those observed in M82 but with 
a much higher luminosity. We show that the starburst region can be divided into an outer zone, beyond a 
radius of about 2.1 kpc, with a gas temperature of about 10 7 K and a central region with temperatures in the 
range (2-6) x 10 7 K. The gas in the outer region emits most of the observed O VIII L a line and the gas in the 
inner region the emission lines of higher ionization ions, including a strong Fe XXV line. We also identify a 
small inner part, very close to the active nuclei, with typical Seyfert 2 properties including a large amount of 
photoionized gas producing a strong Fe Ka 6.4 keV line. The combined abundance, temperature and emission 
measure analysis indicates super solar Ne/O, Mg/O, Si/O, S/O and possibly also Fe/O. The analysis suggests 
densities in the range of (0.07-0.28)£~ 1,/2 cm~ 3 and a total thermal gas mass of ~ 4 x 10 8 e'/ 2 M©, where e is 
the volume filling factor. We used a simple model to argue that a massive starburst with an age of ~ 2 x 10 7 
years can explain most of the observed properties of the source. NGC 6240 is perhaps the clearest case of an 
X-ray bright luminous AGN, in a merger, whose soft X-ray spectrum is dominated by a powerful starburst. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - quasars: emission lines - galaxies: 
starburst 



1. INTRODUCTION 

Ultra-luminous infrared galaxies (ULIRGs) have most 
of their luminosity emerging in the infrared with L(8- 
1000 urn) > 3 x 10 45 ergs" 1 (Sanders & Mirabel 1996). 
Many ULIRGs occur in interacting systems and it is thought 
that galaxy mergers produce concentrations of gas which can 
fuel a giant black hole (BH) and/or result in prodigious star- 
burst activity. ULIRGs are fundamental for the understand- 
ing of the star formation history in the universe, black hole 
merging and growth and as a trigger of activity in galactic nu- 
clei. Yet, the relationship between these systems, "normal" 
active galactic nuclei (AGN) and starburst galaxies is still un- 
clear (e.g. Sanders et al. 1988; Genzel et al., 1998; Levenson, 
Weaver and Heckman 2001a). Of special interest are those 
ULIRGs showing a clear AGN signature in the X-ray but only 
very weak AGN-type emission in the IR. This may indicate al- 
most complete obscuration of the central AGN with important 
implications to the X-ray background and AGN evolution. 

NGC 6240 (z=0.0245; hereafter we assume a luminosity 
distance of 106 Mpc and angular scale such that 1" corre- 
sponds to 490 pc at the source) shows two compact optical- 
IR-radio sources (Tacconi et al. 1999) and two point-like X- 
ray sources (Komossa et al. 2003, hereafter K03). Its IR 
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luminosity, using the Sanders and Mirabel (1996) definition 
and the IRAS fluxes, is L(8-1000 fltn) ~ 7 x 10 11 L© which 
we assume to represent about 90% of the bolometric lumi- 
nosity. While being a bit below the "standard" ULIRG lumi- 
nosity of 10 12 L©, the source has been considered by many 
to be a member of this group. Tecza et al. (2000) show IR 
imaging spectroscopy of the source and discuss the starburst 
in the inner part. More recent IR data are given by Bog- 
danovic et al. (2003) and by the recent high spatial resolu- 
tion near IR imaging of Max et al. (2004) that shows sev- 
eral point sources, possibly young star clusters, in the inner- 
most region. ASCA (Turner et al. 1997; Netzer et al. 1998), 
BeppoSax (Vignati et al. 1999) and Chandra (K03 and Ptak 
et al. 2003) observations show the presence of an absorbed 
Seyfert 2 type continuum with an absorbing line-of-sight col- 
umn of Nh ~ 1.5 x 10 24 cm~ 2 and several large equivalent 
width (EW) lines. The absorption corrected X-ray luminosity 
is L2-io*.fV ~4x 10 43 erg s _1 (obtained from our own analy- 
sis, depending somewhat on the assumed absorbing column, 
see below). Extrapolation of the unobscured X-ray continuum 
to lower energies, assuming a "typical" AGN spectral energy 
distribution (SED), suggests that some 50% of the bolometric 
luminosity is due to this component, yet there are only weak 
high excitation AGN-type IR lines (Genzel et al. 1998, Lutz 
et al., 2003). 

X-ray observations of several well studied type-II AGN 
(e.g. NGC 1068, Mrk3, the Circinus galaxy, see e.g. 
Kinkhabwala et al. 2002) show emission lines from gas which 
is photoionized by the central non-thermal source. The star- 
burst contribution to the X-ray emission in those sources is 
negligible. Different conclusions about the origin of the X- 
ray gas in many sources that show combined Seyfert-starburst 
signs were reached by Levenson et al. (2001a; 2001b) in their 
extensive study of such sources. Unfortunately, all the Lev- 
enson et al. spectroscopic analysis is based on low resolution 
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TABLE 1 

XMM-Newton OBSERVATION LOG FOR NGC6240 



Obs. ID 


Revolution 


Start time 


Duration 




Good time interval (GTI) [ksec] 










[ksec] 


pn 


MOS1 


MOS2 


RGS1 


RGS2 


0101640101 


0144 


2000 Sep 22 01:38:46 


30.1 


12.3 


23.6 


23.6 






0101640601 


0413 


2002 Mar 12 21:37:24 


19.0 


7.3 


15.5 


16.2 


10.1 


9.8 


0147420201 


0597 


2003 Mar 14 18:06:14 


31.6 


5.1 


10.5 


10.5 


14.3 


14.3 


0147420301" 


0599 


2003 Mar 18 21:01:56 


28.1 








2.9 


2.9 


0147420401 


0673 


2003 Aug 13 10:29:32 


14.1 


7.7 


11.8 


11.8 


8.5 


8.5 


0147420501 


0677 


2003 Aug 21 10:25:09 


31.2 


3.6 






12.0 


12.0 


0147420601 


0681 


2003 Aug 29 11:25:15 


9.2 


1.4 






2.9 


2.9 






Total GTI (ks) 


163.3 


37.3 


61.4 


62.0 


50.6 


50.3 



a EPIC-pn data suffered from telemetry problems, and could not be recovered. 



CCD spectra, mostly by ASCA and ROSAT, which can be 
very problematic. The clearest example is NGC 1068 where 
Kinkhabwala et al. (2002) followed suggestions by Marshall 
et al. (1993) and Netzer and Turner (1997) and showed, very 
convincingly, that the so-called "thermal X-ray gas" is in fact 
photoionized gas in the source. There are also several well 
documented cases where abundances as small as 0.1 solar 
have been derived from highly significant fits to CCD spec- 
tra (e.g. Dahlem et al., 1998) while later works showed that 
the reason was the lack of spectral resolution and not the un- 
usual abundances. Apparently, the improved spectral capabil- 
ity provided by XMM-Newton and Chandra can completely 
change the conclusions regarding the origin of the emission 
line gas in such sources. 

Based on its hard X-ray spectrum, NGC 6240 can definitely 
be classified as a very luminous type-II AGN. Yet Boiler et 
al. (2003) and Ptak et al. (2003) argued for a more com- 
plex model involving a combination of an AGN and a star- 
burst. Given the need for a high resolution X-ray spec- 
troscopy, we have obtained data that enabled us to resolve the 
AGN-starburst confusion in this source. §2 below discusses 
new XMM-Newton observations of NGC 6240 and §3 gives 
detailed spectroscopic analysis and physical interpretation of 
the spatially resolved emission in this source. Finally, in §4, 
we discuss implications to the local starburst as well as to the 
general question of starburst activity in giant mergers. 

2. OBSERVATIONS AND REDUCTION 

Table [l] gives a summary of the XMM-Newton data used in 
this paper. It includes all the relevant observations including 
those published by Boiler et al. (2003, first entry in the ta- 
ble). We reduced all data using the Science Analysis Software 
(SAS v5.3.0) in the standard processing chains as described in 
the data analysis threads and the ABC Guide to XMM-Newton 
Data Analysis. 8 

Unfortunately, all the 2003 observations had very restricted 
visibility windows and had to be split into several short point- 
ings. Also, the observations were taken at times when the 
spacecraft was close to the Earth's radiation belts resulting in 
high background count rates. Several of the observations were 
thus terminated before their completion. Tabled lists the to- 
tal exposure time of each XMM-Newton observation and the 
good time interval (GTI) for each of the instruments. 

The EPIC-pn (hereafter pn) and the two EPIC-MOS (here- 
after MOS) detectors were operated during our 2003 obser- 
vations in the large window mode with the thin filter. The 

8 http://heasarc.gsfc.nasa.gov/docs/xmm/abc 



pn and MOS data were processed using SAS in the standard 
way. Source events were extracted from circular regions of 
radius 30" and background events from an annulus with a ra- 
dius of 30"-87". Epochs of high background event were ex- 
cluded. The known PSF properties suggest that some 10% 
of the flux may be emitted outside of the chosen extraction 
window. We have therefore re-extracted our longest GTI ex- 
posure with twice the above radius and a background region 
far removed from the source. The comparison with the results 
of the smaller extraction region show an insignificant differ- 
ence thus the chosen extraction window has little influence on 
our final results. 

Six pn and four MOS observations, with a total GTI of 37 
and 61 ks, respectively, resulted in useful data. After reduc- 
ing a spectrum from each observation and verifying there is 
no time variability (see below) we combined all spectra from 
each detector into one spectrum. This was done using the 
FTOOL task MATHPHA. We also combined the matching 
rmf and arf files using the tools ADDRMF and ADDARF 
(weighting individual exposures by their respective integra- 
tion times). We then binned the combined spectrum of each 
detector to have at least 25 counts per bin. All data presented 
here are corrected for galactic absorption with a column den- 
sity of 5.9 x 10 20 cirr 2 (Dickey & Lockman 1990) 

The RGS 1 and RGS2 were operated in the standard spec- 
troscopy mode and produced 50 ks of GTI. Background sub- 
traction was performed with the SAS using regions adjacent 
to those containing the source in the spatial and spectral do- 
mains (using the defaults given by the SAS rgsproc tool). The 
spectra were extracted into uniform bins of ~ 0.04 A (roughly 
half the RGS resolution and four times the default bin width) 
in order to increase the signal-to-noise (S/N) ratio. Flux cal- 
ibration was obtained by dividing the counts by the exposure 
time and the effective area at each wavelength. Each flux- 
calibrated spectrum was also corrected for Galactic absorp- 
tion and the two RGS spectra of each observation were com- 
bined into a weighted mean. At wavelengths where the RGS2 
bins did not match exactly the wavelength of the RGS 1 bins, 
we interpolated the RGS2 data to enable the averaging. Af- 
ter verifying that all spectra from different observations are 
consistent with each other, we combined the spectra from the 
six observations listed in Table^into an error-weighted mean 
which is presented and analyzed below. At A > 20 A the RGS 
spectrum is so noisy that no meaningful line intensity can be 
measured and even upper limits on emission line fluxes are 
highly questionable. Those long wavelength lines definitely 
require longer integration. 

The spatial resolution of the XMM-Newton instruments was 
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discussed, in great detail, by Boiler et al. (2003). These au- 
thors show that some 90% of the flux is encircled within a 
radius of about 30" corresponding to about 15 kpc. This is 
larger than the entire X-ray bright image observed by Chan- 
dra (see below) hence all our XMM-Newton related results 
refer to the total flux of the source. 

Our work makes use also of the Chandra data published 
by K03. This observation was obtained on 2001 July 29 us- 
ing the ACIS-S3 CCD (Obs. ID. 1590) and resulted with 
a GTI of 37 ks. We obtained the data from the Chandra 
archive and reduced and analyzed them using CIA03.2.1 
(with CALDB3.0.0) and its standard threads for extraction 
of extended and point like sources. We used various spatial 
and energy filters to extract images and spectra as detailed in 
the following sections. We note that 1" corresponds to 0.49 
kpc at the source and the sizes of all the spatially resolved 
images discussed below are much larger than the instrument 
point spread function (PSF) which contains 80% of the power 
within 0.685". Thus none of the linear sizes and volumes de- 
rived below is affected by more than about 10%; a small un- 
certainty considering the other unknowns discussed below. 

Since our observations were split into various segments that 
were separated by weeks to months, we were able to study 
the time dependent of the source. We have carried out a stan- 
dard analysis using the high energy (less obscured) E > 6 ke V 
continuum. We found no variability to within about 10%. We 
note that Ptak et al. (2003) show the historical (1994-2002) 
light curve of the source where they detect variations by up to 
40%. 

3. SPECTROSCOPIC ANALYSIS 

3.1. A comparison of the spectra of NGC 6240 and M82 

M82 is a classical luminous starburst galaxy that has been 
studied by all X-ray satellites (Read & Stevens 2002 and ref- 
erences therein). The galaxy is known to have a prominent 
superwind most clearly seen in H a images as well as in X- 
ray lines (e.g. Stevens, Read and Bravo-Guerrero, 2003). The 
XMM-Newton based work of Read and Stevens and Stevens et 
al. claim a two region superwind with indications for a hotter 
central region. They also suggested abundances that deviate 
from solar and some contribution to the observed high energy 
continuum from a known point source (possibly a medium- 
size black hole) and perhaps several unresolved X-ray bina- 
ries. The starburst origin of most of the observed X-ray emis- 
sion is clearly indicated by the many iron charge states and 
by the line ratios in several He-like ions (Read and Stevens 
2002). 

Fig.[T]shows the RGS spectra of NGC 6240 and M82. The 
M82 spectrum was obtained from the XMM-Newton archive 
and is based on the data presented in Read and Stevens (2002) 
and in Stevens et al. (2003; 30 ks integration, 25 ks GTI). It 
was extracted in the same way as the NGC 6240 spectra. The 
spectra shown here are normalized to allow easy comparison. 
The similarity of the 5.5-20 A spectrum of the two galaxies 
is striking. All strong emission lines seen in the spectrum of 
M82 are reproduced in the (noisier) spectrum of NGC 6240 or 
are consistent with suspected emission features. As explained 
above, the A > 20 A RGS spectrum of NGC 6240 cannot be 
used for obtaining line fluxes or even meaningful upper limits. 

Line intensities for the two galaxies were measured from 
their RGS spectra by fitting simple Gaussians to strong ob- 
served features. Theoretical wavelengths and zero intrinsic 
widths were used as constraints on the measurements. The 
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wavelength[A] 

FIG . 1 . — Comparison of the XMM-Newton spectra of M82 and NGC 6240. 
The RGS spectra are shown in black and the pn spectra in red. The flux scale 
of the two has been arbitrarily shifted to allow easier comparison. 

results are listed in Table 2 (note that the Read and Stevens 
2002 and the Stevens et al. 2003 papers do not list the line 
intensities used in their analysis of M82). Obviously, the un- 
certainties on the line intensities in NGC 6240 are larger due 
to the worse S/N of its RGS spectrum. An important issue is 
the relative intensity of the various line components (r, i and 
/) in He-like ions. These are very good temperature indicators 
provided the effect of line pumping (sometimes referred to as 
continuum fluorescence) is taken into account. Unfortunately, 
the RGS resolution is insufficient to provide a clear separation 
of the various He-like line components for ions with Z > 12. 
Moreover, the S/N of our NGC 6240 RGS spectrum requires 
heavy binning thus we cannot trust individual line compo- 
nents. The measurements of the Ne IX complex suggests it 
is dominated by the r line, the signature of a high tempera- 
ture gas, but the uncertainty is large and we prefer using other 
types of analysis explained below. 

As evident from Table 2, the intensity ratios of most emis- 
sion lines in the two spectra are very similar. There is some 
tendency for the H-like/He-like line ratio of some elements to 
be somewhat larger in NGC 6240. This is most noticeable for 
the silicon lines. However, the silicon lines region is at the end 
of the RGS response and the spectrum in this region is partic- 
ularly noisy. Thus the uncertainty on the flux of the Si XIV 
line is very large. To improve the line measurement, and in 
order to add lines outside of the RGS range, we have included 
in the comparison the MOS and pn spectra of the two galax- 
ies. We fitted models to the data over the 1.7-10 keV range 
using XSPEC version 11.3.0. The NGC 6240 model consists 
of the following components: a. An absorbed power-law. b. 
A large column density line-of-sight neutral absorber, c. An 
additional power-law which may be either the reflected high 
energy continuum (in which case the two were forced to have 
the same slope) or due to local X-ray sources. The most likely 
contributers in this case are X-ray binaries which are known 
to have harder X-ray continua. The power-law in this case 
was forced to have a flatter slope. An additional possibility 
is a very high temperature thermal continuum mimicking a 
power-law. As shown below, such temperatures are inconsis- 
tent with the temperatures derived from the emission lines, 
d. The following emission lines (all Gaussians with zero- 
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TABLE 2 
Measured line fluxes a 



Line 


vv dVGidigiiii/i 




M82 


Tn c tn l m 
ma Li mucin 


Kj3 b 


1.75 


15.5 ± 1.1 




pn+MOS 


Fe XXV 


1.85 


81.4±7.9 


174 Q+29.9 


pn+MOS/pn 


IV u 


1 Q4 


112 1 ± 7 8 


59.2_ 143 




S xvi c 


4.73 


6.0± —3.0 


pn+MOS/pn 


Sxv 


5.04,5.10 


20.5 ±2.8 


240.8i*j;| 


pn+MOS/pn 


Si xiv 


6.18 


16.5±2.1 


126.0±63.5 


pn+MOS/RGS 


Si xm 


6.65 


26.0 ±1.9 


519.0± 150.0 


pn+MOS/RGS 


Mgxn 


8.42 


20.1 ±8.1 


285.0 ±94.0 


RGS 


MgXI 


9.17 


18.5 ±17.5 


335.0±96.9 


RGS 


NeX 


12.13 


25.4±12.5 


115.1 ±40.3 


RGS 


FeXXI 


12.28 


13.7± 10.7 


100.3 ±30.8 


RGS 


NelX 


13.45 


12.2±4.0 


236.0 ±92.1 


RGS 


FeXX, FeXIX 


13.77, 13.79 




65.2±23.3 


RGS 


FeXVIII 


14.21 




63.7 ±16.4 


RGS 


FeXVII, FeXIX 


15.01, 15.08 


5.6±2.6 


151.0±30.4 


RGS 


FeXVII 


15.26 


9.3±3.1 


64.2 ±22.0 


RGS 


FeXVIII, FeXIX 


16.00, 16.07, 16.11 




63.4±21.3 


RGS 


FeXVII 


16.78 




41.4±15.0 


RGS 


FeXVII 


17.05, 17.10 


5.5 ±2.8 


8 1.0 ±26.0 


RGS 


O VIII 


18.97 


12.1±3.6 


69.9 ±27.3 


RGS 



a In units of 10 15 ergs cm 2 s 1 . 

b Line intensity constrained to 0. 1 14 the intensity of Ka. 

c Flux and uncertainty estimated from integrating above a power-law fit to the local continuum. 
d Left is for NGC 6240 and right for M 82. 



widths which gives the instrument resolution when used in 
XSPEC): 1.865 keV (Sixm), 2.006 keV (Si xiv), 2.461 keV 
(S xv), 2.622 keV (S xvi), 6.400 keV (Fe i-xvii, hereafter 
Ka), 6.700 keV (Fexxv), 6.96 keV (Fexxvi) and 7.100 
keV (Fe I-XVII, hereafter Kj3). All line energies were fixed 
at those rest energies and the Ka/K/3 ratio was forced to its 
known theoretical value. % 2 minimization is obtained by fit- 
ting, simultaneously, the pn and the two MOS spectra. 

The results of the fitting process are as follows: The two 
powerlaws with forced identical slopes gave Y = 1. 70 ±0.01. 
The two independent power-laws gave Y = 1.76 ±0.02 for 
the absorbed one and Y = 1.19 ±0.02 for the ("binary") 
unabsorbed power-law. The absorbing columns are Nh = 
1.37 ±0.02 x 10 24 cm" 2 for the first case and N H = 1.35 ± 
0.03 x 10 24 cm~ 2 for the second case. The slopes and the col- 
umn densities are in rough agreement with the Netzer et al. 
(1998), Vignati et al. (1999) and Ptak et al. (2003) earlier re- 
sults. In all cases, acceptable fits were obtained (with reduced 
X 2 in the range 0.79-0.92). The observed EWs for the high- 
est energy lines are 1.09 keV for the Ka line and 0.81 keV 
for the Fe XXV line. Both fits do not require any Fe XXVI line 
although we cannot exclude a weak line that contributes up to 
30% to the feature identified here as K/3 . 

We have made no attempt to fit the entire 0.3-10 keV pn 
and MOS spectra of NGC 6240 in a way similar to several 
earlier works. Such models are necessarily limited by the low 
spectral resolution of CCD-type detectors that results in a de- 
generacy with respect to the composition and temperature of 
the emitting gas. As explained in §1, several earlier works us- 
ing such fits resulted in erroneous conclusions. Our findings 
are based as much as possible on the RGS data at long wave- 
lengths, where such a degeneracy does not exist. The pn-MOS 
fit is therefore only at high energies. Nevertheless, we believe 
that the Si XIII and Si XIV line fluxes obtained from such fits 
are more reliable than the ones obtained from the RGS spec- 



.0012 rn 



Fe XXV Ka 

.001 — 



^.0008 — 
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FIG. 2. — The iron-K complex in NGC 6240 (combined pn and MOS spec- 
tra). 



trum, which is very noisy at those energies. Therefore they 
are the ones listed in Table 2. 

A small part of the pn spectrum obtained with the tabulated 
line intensities is added to Fig.^on the same scale as the RGS 
spectrum. The iron-K region (unfolded spectrum, combined 
pn and MOS data with the line fluxes listed in Table 2) is 
shown in Fig. [2] 

The fit to the pn spectrum of M82 was done for the pur- 
pose of measuring lines with wavelengths shorter than 6 A. 
This was achieved by fitting zero-width lines to a continuum 
that was determined in two different ways: a. a single pow- 
erlaw. b. the multi-temperature model of Read and Stevens 
(2002) with zero abundance for all heavy elements. The re- 
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suits of the two methods are similar and the line intensities 
reported in Table 2 are the ones obtained with the second con- 
tinuum. All other lines, with longer wavelengths, were mea- 
sured by fitting local continua and Gaussian-shaped lines to 
the RGS spectrum of the source as done for NGC 6240. We 
note that the purpose of our work is to model the spectrum of 
NGC 6240, hence the only aim of discussing the M82 spec- 
trum is to provide a basis for comparison for NGC 6240. 

3.2. Thermal gas in NGC 6240 

The great similarity of the 4-20 A spectra of NGC 6240 
and M82 suggests that the large scale X-ray emission in 
NGC 6240 is dominated by a massive starburst. Given the 
strong Neix and FeXVll lines, that are typical of low temper- 
ature plasma, and the prominent silicon and sulphur H-like 
lines typical of much higher temperature gases, we suggest 
that there are at least two regions with very different temper- 
atures in this source. Boiler et al. (2003) reached a simi- 
lar conclusion by comparing their XMM-Newton observations 
with the spectrum of NGC 253, another well known starburst 
galaxy. These findings are also supported by the multi-energy 
maps of K03 that suggests that the gas emitting the 0.5-1.5 
keV flux is mostly in the outer parts of the system, further 
away than 2-3 kpc from the two nuclei. According to them, 
the gas emitting the 1.5-2.5 keV radiation is closer in and 
there is a small inner part that emits most of the E > 5 keV 
radiation. 

The parameters of the models presented below are the tem- 
peratures, the abundances of oxygen, neon, magnesium, sili- 
con, sulphur and iron, and the emission measure, 

EM{T) = J n e n H edV cm" 3 , (1) 

where e is a volume filling factor. We have calculated sev- 
eral low density thermal plasma models for a range of tem- 
peratures and compositions and compared them with the line 
list given in Table 2. All models shown here, as well as 
the photoionization models discussed in §3.3, were computed 
with ION2004, the 2004 version of the code ION (Netzer 
et al. 2003 and references therein). The definition of "so- 
lar composition" as used in these models for the elements 
under study here is N(0)/N(H)=5 x 10~ 4 , N(Ne)/N(H)=l x 
10~ 4 , N(Mg)/N(H)=3.5 x 10~ 5 , N(Si)/N(H)=3.5 x 10~ 5 , 
N(S)/N(H)=1.6 x 10~ 5 andN(Fe)/N(H)=4 x 10~ 5 . An impor- 
tant addition to the earlier version of the code is the inclusion 
of more realistic collision strengths for various iron ions and 
the better treatment of the various processes in FeXVll (Doron 
and Behar 2002). 

The first stage of the analysis includes re-examination of 
the spatial distribution of the emitting gas. This was done 
by separating images obtained from the Chandra archive into 
narrow energy bands around strong lines and by extracting 
Chandra ACIS-S spectra. We found that most of the flux in 
the O VIII line is emitted in the outskirt of the source, in a 
few (4-5) localized regions outside a radius of about 2. 1 kpc. 
These "lobes" are clearly seen in several of the K03 images. 
The Ne X line originates, in almost equal amounts, in the lobes 
and in a spherically looking central part with a radius of 2.1 
kpc. We identify one component of the thermal gas model 
with the lobes and name it "the outer zone". The temperature 
of this gas is probably between the peak emissivity temper- 
atures of the O VIII and Ne X lines (3 - 8 x 10 6 K). We also 
identify a higher temperature, spherical looking zone showing 
some Ne X emission and most of the Mg XII, Si XIV and S XV 
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FIG. 3. — Thermal line emission as a function of electron temperature 
for solar composition plasma. The curves show computed L a line fluxes for 
various H-like ions relative to the peak emissivity of the O VIII L a line. Points 
with error bars give the observed ratios in NGC 6240. 

line emission. By combining all pixels showing Ne X emis- 
sion outside of the 2. 1 kpc radius, we can evaluate the fraction 
of the volume emitting this line. Our two dimensional line 
emission maps translate to roughly 115 kpc 3 which is about 
half the volume between 2.1 and 4 kpc. The inner zone that 
shows emission of Ne X and higher ionization species has a 
volume of about 39 kpc 3 . 

Second, we examined the intensity ratios of the strongest 
Fe XVII lines relative to O VIII. The emissivities of these lines 
peak very close to the emissivity of the Ne X lines, but unlike 
the O VIII and Ne X lines, it drops sharply above about 10 7 K. 
We found that at temperatures of order 10 7 K, and with Fe/O 
solar to twice solar, a collisional model gives a satisfactory 
fit to the relative intensities of the Fe XVII and O VIII lines. 
Thus the thermal gas model of NGC 6240 requires between 
solar and 2 x solar Fe/O. We did not attempt any model where 
different temperature components have different metalicity as 
the data do not warrant such an attempt. 

Third, we investigated the higher temperature component 
which is required to explain the strong Si XIV, S XVI and pos- 
sibly Fe XXV lines. The observed line ratios reflect the met- 
alicity, temperature and EM of the gas in the various regions. 
This can be illustrated with the help of the generic solar com- 
position model shown in Fig. [3] The figure shows all the H- 
like L a lines observed by us relative to the peak emission of 
the O VIII L a line as a function of temperature. This is the 
best way to illustrate the intensities of collisionally excited 
lines that are emitted at temperatures close to their peak emis- 
sivity temperature. Obviously, the real situation in NGC 6240 
is more complicated but the diagram helps illustrate the main 
point regarding the need for temperature distribution in the 
source. 

As evident from Fig. [3] the assumption that most line emis- 
sion occurs at temperatures close to the peak emissivity tem- 
perature of the various ions clearly indicates that neon, mag- 
nesium, silicon and sulphur require either larger abundances 
and/or much larger EM in their production regions compared 
with oxygen. This confirms the findings based on the imag- 
ing analysis that there must be an inner central region with 
temperatures that are large enough (T > 1.5 x 10 7 K) to ex- 
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plain the intensities of the magnesium, silicon and sulphur 
lines. If the temperature in some parts of this region exceeds 
~4x 10 7 K, it can also explain the observed Fe XXV line. 

As explained, the line emissivities shown in Fig. [5] suggest 
possible over abundances of Ne/O, Mg/O, Si/O and S/O. Be- 
fore investigating this in detail we note that there are indica- 
tions for high abundances in M82, our comparison source (the 
abundances found by Read and Stevens (2002) for M82 are 
not entirely consistent with the Stevens et al. (2003) results 
although the two investigations are based on the same XMM- 
Newton observation). This is also the conclusion of the recent 
Origlia et al. (2004) analysis of the source. We have no way 
to scale the metal abundances in NGC 6240 relative to hydro- 
gen because of the very large uncertainty on the intensity of 
the low temperature thermal continuum. Thus, all the abun- 
dance analysis presented here is relative to oxygen. Having 
no way to deduce the oxygen abundance, we assumed solar 
O/H. While none of our conclusions depend on this assump- 
tion, we note that it is consistent with the Read and Stevens 
(2002) finding for O/H in M82. 

We carried out a combined abundance and temperature 
analysis aimed to find the combination of temperature, EM 
and composition that best fits the spectrum of NGC 6240. 
Given the large uncertainty on most line fluxes, we only in- 
vestigated three zone models and four different compositions 
ranging from solar to 2.5 x solar metalicity. We treat the el- 
ements in two groups, one containing Ne, Mg, Si and S, 
and one containing Fe. We have constructed a grid of ther- 
mal models covering the temperature range 10 6 - 10 8 K and 
changed the assume EMs for the various zones over a large 
range. The results were tested against seven observed line ra- 
tios (seven line intensities relative to the O VIII 18.96A line): 
L a of H-like neon, magnesium, silicon, and sulphur, the He- 
like Fe XXV line, the Fe XVII line at 15.0 A and the upper limit 
on Fe XXVI. The results of the analysis are: 

1 . No solar abundance model can produce all the observed 
line ratios with their uncertainties. The best solar com- 
position model can fit only 4 of the seven line ratios. 

2. Models with 1 .5 - 2.5 x solar for Ne/O, Mg/O, Si/O and 
S/O can fit all line ratios to withing the uncertainties 
provided Fe/O is between solar and 2.5 x solar. The 
higher abundance cases (2 and 2.5 x solar) produce sig- 
nificantly better fits and a larger variety of acceptable 
combinations of temperatures and EMs. 

3. The spatial distribution of the Ovill emission can be 
used as an additional constraint on the model. The as- 
sumption that at least 60% of the total line emission 
comes from the outer zone (as our analysis clearly indi- 
cates) does not change any of the previous conclusions. 
Assuming this fraction is more than 85%, which is also 
consistent with the available data, forces the tempera- 
ture of the outer region to above 10 7 K and gives poorer 
fits mostly for the Fe XVII lines. Unfortunately, accu- 
rate determination of this fraction using the available 
data is difficult to obtain. For example, we cannot ex- 
clude the possibility that some of the O VIII emission 
originates even further from the center than indicated 
by the Chandra images. Such a component is likely to 
have temperatures much below 10 7 K. 

4. There is no unique solution to the temperature-EM dis- 
tribution in the source. In fact, the temperatures and 
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FIG. 4. — A three-zone model fit to the spectrum of NGC 6240. The black 
and red lines are the RGS and pn spectra shown in Fig. 1 and the blue line is 
the three-temperature model. The RGS data is binned to 0.08A and the model 
assumes Gaussian emission lines with FWHM corresponding to to this width. 

EMs in the three zones can change by up to 0.1 dex 
while still being inside the range allowed by the ob- 
servations. Thus, there are dozens of solutions within 
these uncertainties and no statistically sound way to tell 
them apart. 

The conclusion of the analysis is that a multi-temperature, 
multi-zone thermal plasma model can explain all observed 
emission lines except for the Ka and Kj3 lines, to a reasonable 
accuracy. Table [3] shows two sets of EMs and temperatures 
that successfully reproduce all seven line ratios for a case of 
solar Fe/O and twice solar for the other elements, and a sec- 
ond case of 2 x solar for Ne/O, Mg/O, Si/O, S/O and Fe/O. 
Table|4]shows the calculated line ratios for the two cases and 
compares them with the observations. In the rest of the paper 
we focus on the first case (model A in the table). Fig.|4]shows 
a fit of this model to the 4-20 A spectrum of NGC 6240. 

There are two fundamental uncertainties related to the sug- 
gestion that the 2.1 kpc radius region can be separated into 
two sub-zones with the temperatures found earlier. They are 
both related to the role of photoionization close to the two ac- 
tive nuclei. The first is the origin of the FeXXV 1.85 A (6.7 
keV) line. This line can be produced by hot thermal plasma, 
as shown in our calculation, or by recombination in a colder 
photoionized gas, as seen in several well studied type-II AGN. 
The latter case is investigated in §3.3. The second uncertainty 
is related to the relatively strong 2-5 keV power-law contin- 
uum that is required by the spectral fit. As explained in §3.1, 
such a continuum can be produced by a very high temperature 
gas, by a local population of X-ray binaries or by reflection of 
the central power-law continuum. The second possibility has 
been discussed in several recent papers, see e.g. Persic and 
Rephaeli (2002). The predicted continuum in this case is very 
hard, with r ~ 1.2. The third mechanism is known to domi- 
nate the soft X-ray emission of other type-II AGN at energies 
where the central source is obscured. In such cases, the re- 
flected continuum is seen from an extremely small, usually 
unresolved part of the source. 

To study this issue we have extracted the Chandra spectra 
of three separate regions of the source: the outer zone be- 
tween 2. 1 and 4 kpc, the central zone with a radius of 2. 1 kpc, 
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TABLE 3 

Temperatures (in K) and emission measures (in cm~ 3 ) for thermal models 



Case 


Ti 


EM\ 


T 2 


EM 2 


T 3 


EM 3 


Model A 
Model B 


10 7 
8x 10 6 


1.9 x 10 64 
1.2 x 10 64 


1.6 x 10 7 
1.6 x 10 7 


1.4 x 10 64 
2.0 x 10 64 


6.3 x 10 7 
5 x 10 7 


5.4 x 10 64 

3.5 x 10 64 



TABLE 4 

Line ratios relative to O viii X 18.966 for the models in Table0 





O VIII 


Nex 


MgXII 


Si XIV 


Sxvi 


FeXXV 


Observed 


1.0 


2.10± 1.21 


1.66 ±0.83 


1.36 ±0.44 


0.50 ±0.30 


6.73±2.10 


Model A 


1.0 


1.60 


1.14 


1.22 


0.56 


5.99" 


Model B 


1.0 


1.55 


0.93 


1.16 


0.53 


6.26 b 



"Solar Fe/O and twice solar for Ne/O, Mg/O, Si/O and S/O. 
"Twice solar for Ne/O, Mg/O, Si/O, S/O and Fe/O. 
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FIG. 5. — Chandra spectra of the outer (blue), central (red, all emission 
within 2.1 kpc from the center) and innermost (black) zones in NGC 6240. 
Note that most of the 4-6 keV continuum comes from the innermost part 
while the Fe XXV line seems to originate mostly in the central region. 

and an innermost zone encompassing the two nuclei. The lat- 
ter is seen as an elliptical region and can be described by an 
ellipsoid with a major axis of 1 .4 kpc and two other axes of 
0.7 kpc each. The larger dimension is along the line connect- 
ing the two nuclei and the separation between the nuclei is 
roughly the size of the semi-major axis. The corresponding 
volume is about 0.36 kpc 3 . The high energy parts of these 
spectra are shown in Fig. [5] 

The comparison of the three spectra confirms the general 
suggestion that the outer zone contributes most to the low 
energy emission while the inner zones produce spectra that 
are typical of a much higher temperature gas. Unfortunately, 
some features in the Chandra spectra are different from those 
seen in the combined pn+MOS spectra for reasons that we do 
not understand. This is seen mostly as excess emission be- 
tween about 2.1 and 2.4 keV. The problem has already been 
noted by Ptak et al. (2003) in their modeling of the K03 
data and was suspected to be a calibration error due to a well 



known mirror feature near the gold edge. We are not sure how 
to confirm this idea since we used the best available calibra- 
tion provided by the CXC (see §2). Given this uncertainty, 
we decided not to use the Chandra spectra for measuring line 
intensities and only to address, qualitatively, the issues of the 
Fe XXV line and the reflected continuum. 

As seen in the full Chandra spectra (not shown here), and 
as hinted in Fig.|5J there is a tendency for the relative intensity 
of the 2-5 keV continuum to increase inward, supporting the 
idea of reflection in the photoionization region. The opposite 
is true for the Fe XXV whose flux seems to originate in an in- 
termediate part outside the elliptical photoionization region. 
These findings are marginal because of the difficulty in sep- 
arating the two small zones in the center and because of the 
problematic Chandra spectra. They are shown here in con- 
junction with the theoretical argument presented below (§3.3) 
and seem to support it. The situation regarding the Ka line is 
clearer and the diagram shows that the line originates in the 
inner zone. The small elliptical zone is termed below "the 
photoionized region". 

Given the qualitative impression based on the Chandra 
imaging and Fig. [5] and the detailed modeling of the ther- 
mal gas, we suggest the following three zone structure for the 
starburst dominated regions in NGC 6240. The first part is the 
one addressed earlier, outside a 2.1 kpc radius. The electron 
temperature there is 8 — 10xl0 6 K and the volume occupied 
by the line emitting gas is about 115 kpc 3 . Inside the 2.1 
kpc radius there is an inner spherical part (hereafter "the in- 
ner zone") where the temperature is 5 — 6 x 10 7 K and where 
most of the Fe XXV is emitted. Outside of this zone there is an 
outer shell (hereafter the "intermediate zone") where the tem- 
perature is ~ 2 x 10 7 K and where most of the magnesium 
and silicon line emission takes place. We also suggest that 
the small-volume innermost photoionized zone is part of the 
inner thermal plasma zone. We adopt, somewhat arbitrarily, a 
radius that divides the 2.1 kpc sphere into two concentric re- 
gions of equal volume (19.5 kpc 3 ). This division is not based 
on emission line imaging and is, therefore, rather uncertain. 
However, the general abundances derived here, and the total 
hot gas mass deduced below, are rather insensitive to the exact 
division. 
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FIG. 6. — Line ratios relative to the O VIII L a line for photoionized gas 
with twice solar abundance for Ne/O, Mg/O, Si/O and S/O and solar O/H 
and Fe/O. The column density is 10 22 cm~ 2 except for the Fe-Ka line where 
the column density is 10 23 cm~ 2 . Points with error bars, or upper limits, 
represent the data in Table 2 and are arbitrarily located along the horizontal 
axis to illustrate the difficulty to obtain the observed line ratios for a specific 
value of U (oxygen) . 

3.3. Photoionized gas in NGC 6240 

The two strong point-like X-ray continuum sources ob- 
served in NGC 6240 can ionize low density gas many pc away 
from the center. Such gas will produce fluorescence and re- 
combination iron lines. The intensity of the lines depends on 
the covering factor and the column density of the illuminated 
gas and its level of ionization. The strong 6.4 keV Ka iron 
line (Fe I-XVII) observed in NGC 6240 cannot be produced 
by the starburst gas and must be due to fluorescence. Indeed, 
the K03 findings and our own work (Fig0 are consistent with 
the production of this line only in the innermost photoionized 
part of the source. In contrast, the FeXXV 1.85 A line can 
be produced either by hot gas in starburst regions or in high 
ionization photoionized gas. A good example of the later is 
the case of NGC 1068 (Marshall et al., 1993; Netzer & Turner 
1997; Kinkhabwala et al., 2002). Below we discuss this pos- 
sibility for the special case of NGC 6240. 

Fig-EJshows the results of photoionization calculations that 
predict the intensities of the Nex, MgXll, Si XIV, Sxvi, 
Ka, Fe XXV and Fe XXVI lines relative to the intensity of the 
O VIII L a line. The intensity ratios are shown as a function of 
the ionization parameter U (oxygen) which is the photon to the 
hydrogen density ratio defined over the energy range 0.54-10 
keV (e.g. Netzer et al. 1998). The calculations were done for 
the composition found here (model A), a power-law contin- 
uum with r(0. 1—50 keV) = 1 .8, a small UV bump (the exact 
shape of which is not relevant to the present calculations), and 
a range of column densities between 10 20 and 10 23 cm~ 2 . The 
diagram shows the 10 22 cm~ 2 case but the main conclusions 
are the same for all. 

The new calculations show that many line ratios involving 
H-like ions are in rough agreement with the observations of 
NGC 6240 given enough freedom in the choice of U (oxygen). 
This is not surprising since the most important temperature di- 
agnostics for the elements considered here, the r and / lines 
of the He-like ions, are missing from the diagram. The f/r 
ratio can, by itself, separate thermal from photoionized gas. 



As explained earlier, the combination of poor resolution and 
S/N prevented us from using such ratios in a meaningful way 
thus we resort to alternative methods. The next important con- 
clusion is that the calculated Fe XXV line (the sum of all line 
components) is too weak relative to its observed value even if 
we allow for a higher iron abundance. Furthermore, the calcu- 
lations show that photoionized gas producing strong Fe XXV 
must also produce strong Fe XXVI which is in contradiction 
with the observations (see Fig.|6jl. Thus, the photoionization 
origin of the Fe XXV line in NGC 6240 is practically ruled out 
by the calculations. 

We also considered the production of Ka and K/3 lines. 
These were computed in a similar fashion and the Ka line 
is shown in Fig. [6] except that in this case we give its rela- 
tive intensity for a column density of 10 23 cm~ 2 (the smaller 
columns cannot produce the observed line intensity). Obvi- 
ously, the line is produced in a much lower ionization gas 
which is entirely consistent in terms of flux and EW with 
dozens of type-II observations and with theoretical calcula- 
tions (Netzer et al. 1998 and references therein). Thus, the 
relative intensity of the lowest and the highest ionization iron 
lines, plus the success of the spatially resolved thermal gas 
model, all suggest that in NGC 6240, Ka and Kj3 are the only 
strong emission lines produced in the photoionized gas. 

4. DISCUSSION 

The new data and modeling of NGC 6240 presented in this 
paper show an unusual ULIRG which is perhaps the clearest 
case, so far, of a very luminous type-II AGN combined with 
a very high luminosity X-ray starburst. The outer parts of the 
source, beyond about 2.1 kpc, are clearly dominated by the 
starburst which is very similar in its properties to the starburst 
observed in M82, except for the much higher luminosity. The 
innermost part, very close to the two active nuclei (an ellip- 
soid of about 0.7 x 0.7 x 1.4 kpc), shows clear signatures of 
photoionized gas; a strong Ka 6.4 keV line and possibly a re- 
flected continuum. The region between the photoionized part 
and a radius of about 2. 1 kpc is dominated by a high tempera- 
ture starburst with strong emission lines of H-like and He-like 
magnesium, silicon and sulphur and He-like lines of Fe XXV. 
The analysis suggests a temperature gradient towards the cen- 
ter which we model as a three zone starburst region. 

NGC 6240 is among the brightest known X-ray luminous 
starbursts. The integrated 0.5-5 keV luminosity, plus the lu- 
minosity of the Fe XXV line, is about 1 .3 x 10 42 erg s _1 . This 
is about a factor 50 more luminous than the starburst in M82. 
There are several other ULIRGs with comparable soft X-ray 
luminosities. For example, Franceschini et al. (2003) list two 
ULIRGs with 0.5-2 keV luminosity exceeding 2 x 10 42 ergs 
s _1 . However, according to the same authors, most of the X- 
ray luminosity in those objects is due to a power-law source 
which is likely to be X-ray emission from numerous binaries. 
The thermal gas emission in those cases is less than we find 
for NGC 6240 (see their Table 3). It is obvious that high reso- 
lution X-ray spectroscopy is key to the understanding of such 
sources and, so far, our data set for NGC 6240 is the only 
one that allows clear separation of the AGN and the starburst 
components. 

The emission measures found here (Table 3), combined 
with the volumes and temperatures of the three thermal-gas 
zones, can be converted to density and mass of the ionized 
gas within the uncertainties on the filling factor. The latter is 
a problematic issue since there are no direct measurements of 
the filling factor for this or other starburst galaxies. Hydrody- 
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namical simulations by Strickland and Stevens (2000) suggest 
that for X-ray bright galactic winds, 10~ 3 < £ < 10 -1 . How- 
ever, larger values are not excluded by any observation known 
to us. In the following we assume e ~ 0.1 but keep in mind 
the large uncertainty on its value. 

Using the model A parameters from Table we get 
n^e 1 / 2 ~ 0.07 cm~ 3 for the outer zone, n^e 1 / 2 ~ 0.14 cm~ 3 
for the intermediate zone and n H e l l 2 ~ 0.28 cm for the 
innermost zone. Assuming e = 0.1 we get densities in the 
range 0.2-1 cm~ 3 , similar to those estimated by Strickland 
and Stevens (2000) for the warm component in several of their 
wind models. The mass of the hot gas is given by 

M hot oc EMn H x oc e l/2 EM l/2 V 1/2 . (2) 

This gives for model A M hot ~ 1.9 x 10 8 e 1 / 2 M© in the outer 
zone, Mj wt — 6.8 x 10 7 e 1//2 M© in the intermediate zone and 
Mhot — 1.3 X 10 8 £ 1/2 M e in the inner zone. Under these as- 
sumptions, and for a uniform value of e — 0.1 in all zones, 
the total mass of thermal gas is about 1.2 x 10 8 M©. The to- 
tal mass derived assuming model B is about a factor of 1.15 
smaller. 

The mass of hot gas obtained here can be compared with 
the expected properties of a young system (e.g. colliding 
galaxies) undergoing violent star formation. To estimate those 
properties we examined a starburst model computed by the 
code STARS (A. Sternberg, private communication). The 
model assumes a Kroupa (2001) IMF and an exponentially 
decaying star formation rate (SFR) with a characteristic time 
of 20 mega-years (Myr). The SFR is adjusted to give the bolo- 
metric luminosity of NGC 6240 at its peak which occurs after 
about 10 Myr. However, the peak is very broad and the lumi- 
nosity is larger than 80% of the peak luminosity at all times 
between about 4 and 20 Myr. The mass of cold gas converted 
to stars in the first 10 Myr is about 2.5 x 10 9 M© and the gas 
ejected by type-II supernovae is about 3 x 10 7 M© after 6 Myr, 
about 10^ M© after 10 Myr and about 3 x 10 8 M© after 20 
Myr. All properties predicted by the theoretical model for the 
time interval of 10-20 Myr are in general agreement with the 
properties deduced by our thermal plasma analysis assuming 
£ = 0.1. 

The comparison of the observed and predicted hot gas mass 
is only valid if the ejected gas stays in the system for the en- 
tire duration of the burst. This may not be the case for the sys- 
tem under study. Assuming a mean outflow velocity for the 
thermal gas of v out km s and a typical radius of 4 kpc, we 
find a dynamical time of about (4000 km s /v ut) Myr, i.e. 
of the same order as the duration of the burst. NGC 6240 is 
known to have a massive superwind associated with the cen- 
tral starburst with a measured H a outflow velocity of about 
1000 km s~ : (Heckman, Armus and Miley, 1990). This may 
indicate similar velocities for the X-ray gas in the source. In 
such a case, our deduced hot-gas mass is probably a lower 
limit to the total ejecta by the starburst provided most of the 
hot gas mass is carried outward by the wind (a rather uncer- 
tain assumption, see Strickland and Stevens, 2000). While it 
is not our intention to fully model this complex scenario, we 
note that the uncertainty on the filling factor is large enough to 
accommodate this case. Incidentally, the dynamical timescale 
of the two nuclei in NGC 6240 is of the same order, about 7 
Myr (e.g. Tecza et al. 2000). 

The results of the abundance analysis are more difficult to 
interpret. Our model suggests super-solar metalicity relative 



to oxygen. This is problematic for a-process elements and 
may indicate that the model is over-simplified. The global 
properties are, however, in general agreement with the sim- 
ple starburst model we have used. In particular, given type-II 
supernovae explosions as the main source of enrichment, the 
composition of the ejected gas after about 10 Myr must be 
solar to within a factor of 2 for all heavy elements, including 
iron. This is consistent with our general finding. 

Tecza et al. (2000) published the results of a detailed IR 
imaging spectroscopy of NGC 6240. They found that the K- 
band light is dominated by red supergiants formed in the two 
nuclei in starbursts triggered about 20 Myr ago. They fur- 
ther deduced, based on their Brackett— y observation, a very 
short duration burst (~ 5 x 10 6 ys) which is already decaying. 
Despite the somewhat different starburst parameters of their 
model, the general results are in rough agreement with our X- 
ray analysis regarding the age of the burst and the mass of the 
ejecta. The finding of a very powerful starburst also agrees 
with the ISO spectroscopy of Lutz et al. (2003). Pasquali, 
Gallagher and de Grijs (2004) have modeled HST observa- 
tions of NGC 6240 and constructed starburst models for the 
various parts of the source. Their conclusions regarding the 
age and the star-forming rates are similar to those found by 
us. To conclude, our findings are consistent with a massive 
star forming event with an age of 10-20 Myr. It remains to be 
seen whether other X-ray bright ULIRGs can give similar in- 
dications or whether NGC 6240 is an unusual member of the 
group. 

Finally we note two general uncertainties in this and sim- 
ilar sources. The separation of the source into a highly ob- 
scured AGN and completely exposed starburst regions is, of 
course, simplified. Some of the innermost highest temper- 
ature starburst regions can also be obscured which will af- 
fect the observed line intensities. We avoided this issue in the 
present study since using a global absorption correction fac- 
tor is meaningless and a proper analysis requires much better 
mapping of individual starburst regions. This requires deeper 
Chandra imaging of NGC 6240. 

The second uncertainty is also related to the unknown ge- 
ometry of the starburst gas. Our modeling assumes that the 
hot thermal gas is not exposed to the central nonthermal ra- 
diation. This can be justified by the large line-of-sight ob- 
scuring column. However, some geometries would allow the 
starburst gas to see this radiation which can affect the level 
of ionization and the emitted flux in the starburst gas. In par- 
ticular, at E> 5 keV, the cold gas may be transparent in sev- 
eral directions and thus the observed intensity of the Fe XXV 
line can be affected by continuum pumping which was not 
included in our modeling of the highest temperature gas (but 
was included in the photoionization calculations). More com- 
plex models that include all those processes must await better 
Chandra and XMM-Newton observations of the source. 
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